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An Inherent area of interest in high-speed flight such as 
experienced in the re-entry of the space shuttle is the creation 
of unusual aerodynamic fields around the space vehicle and its 
associated problems. High pressure fluctuations in the aero- 
dynamic field due to attached and detached boundary layers and 
shock waves give rise to high structural vibrations as well as 
intensive acoustical atmospheres in the interior of the vehicle. 

The problem is two-folded. First, the creation of the boundary 
layers and shock waves is considered as an input to the struc- 
ture of the space vehicle. Second, with the input, the response 
of the structure can be obtained. The determination of the 
structural response however, is complicated due to the random 
pressure loading from the turbulant boundary layer and by the 
acoustical transmission phenomenon. In other words, as the struc- 
ture vibrates, energy is transmitted to the interior of the space 
vehicle, and the rate of energy transmission and sound intensity 
are governed by the mechano-acoustical boundary condition of the 
interior structural boundary. One simplication is often taken 
here. This is to ignore the coupling effects and solve for the 
structural vibration based on the input. The interior acoustical 
condition is solved separately by considering sound transmission 
through the structure from an acoustical input equivalent to the 
input pressure fluctuation within the boundary layer. 

The determination of the pressure fluctuation within the 
boundary layer has been investigated previously by many researchers. 
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Publications are numerous including theoretical treatments by 
Lighthill (1), Phillips (2), and Ribner (3) and experimental in- 
vestigations by Chyu and Manly (4), Coe (5) and Bull and Willis 
(6). Through these efforts, it can now reasonably assert that 
the input pressure fluctuation is known. 

The solution of the structural vibration is not so assured. 
This is due to the requirement for light weight material with 
great strength and stiffness for the structure. New materials 
are contemplated whose behavior under dynamic loading has not 
been fully understood to warrant their application in actual 
flight. One promising material is the composite materials such 
as epoxy reinforced by carbon fibers. It is interesting to stu- 
dy the effect of those parameters such as modulus of elasticity 
(ratio of fiber to matrix), fiber size and orientation, matrix 
strength, dissipative characteristics, etc., so that optimum 
structural design for vibratory input can be obtained. The same 
material is also examined for noise transmission characteristics. 
The tasks for the present grant can be summarized as follows: 

1. The investigation of elastic wave propagation and attenuation 
in a model fiber - matrix . 

2. The measurement of damping characteristics in graphite epoxy 
composite material. 

3. The construction of a sound transmission test facility suit- 
able to incorporate into NASA Ames wind tunnel for measure- 
ment of transmission loss due to sound generation in boundary 
1 ayers . 

4. The measurement of transmission loss of graphite epoxy com- 
posite panels. 
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RESULTS 


The results of the investigation of v.ie task assignments 
mentioned in the last section are briefly mentioned here. More 
formal and complete results in the form of journal publications 
and other presentations are given in the Appendix. 

1. Elastic Wave Propagation in a Model Fiber-Matrix 

The model chosen for this experiment was a metal plate fitted 
with cylindrical inclusions of a different metal to resemble fiber 
elements in a homogeneous matrix material. The input was simulated 
by a sharp pulse executed by explosive or the impingement of a sharp 
knife-edged ram. The application of the pulse was at the edge of 
the plate such that the ensuing elastic waves propagated in a direc- 
tion perpindicular to its inclusions. The propagation of the waves 
were monitored by strain gages. 

Theoretical analysis was performed using Fermat's principle 
of the ray theory where the location of the wave fronts with res- 
pect to time was obtained. The composite solution was compared 
to experimental results and was found to correlate well. This 
presented concrete evidence that the analytical method can be 
used in fiber reinforced composite material analysis. 

This part of the project has supported one graduate student 
to complete his MSME thesis entitled "Theoretical and Experimental 
Investigation of Stress Wave Attenuation in Model of Fiber Rein- 
forced Composite Material" (by Y. Ueng, 1972) and resulted in an 
AIAA Journal publication (Volume II, No. 4, pp. 472-477, April 1973) 
by its Principal and Co-principal investigators. A copy of the 



4 


reprint of the journal article is included in the Appendix. Pre- 
liminary result of the investigation was also presented in the 13^” 
AIAA/ASME/(SAE Structures, Structural Dynamics and Material Conference, 
San Antonio Texas, April 10-12, 1972). 

2. Measurement of Damping Characteristics in Graphite Epoxy Materials 

A technique called random decrement was developed which makes 
possible the computation of damping values as well as the detection 
of damage in structure when only response data is available. Using 
this technique, damping ratios were computed for several modes of 
randomly excited graphite-epoxy panels. The random input was in- 
duced by speakers within one of the twin reverberation chamber built 
for this project. Response data were obtained from miniature accele- 
rometers attached to the panel . 

This part of the project has supported one graduate student to 
complete his MSME thesis entitled “The Measurement of Damping and 
the Detection of Damage in Structures by the Random Decrement Tech- 
nique" (by D. Caldwell, 1975) and re'^ilted in a paper entitled 
"Measurement of Damping in Composite Beams and Plates." Published 
in the 46^*^ Shock and Vibration Bulletin, Nov. 1975, with D. Cald- 
well. See Appendix. The result of research indicated that the 
damping ratios compared well to those which were derived from 
power spectral density method. Further work on the response of 
plates due to random pressure loading from a turbulent boundary 
layer was also coordinated which resulted in a presentation in the 
Navy-NASTRAN Colloqueum, Bethesda, Maryland, September 1974 and 
published in the Proceedings of the Fifth Navy-NASTRAN Colloqueum, 
pp.81-86. See Appendix. 
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3. Sound Transmission Test Facility 

The sound transmission test facility was designed to attach 
to one of NASA Ames wind tunnel to allow for sound transmission 
loss measurements through a composite panel. One side of the 
panel was designed to fit flush with the interior of the wind 
tunnel and become a part of the walls. Sound generation within 
the boundary layer was transmitted through the composite panel 
and the transmission loss was obtained in sites by this manner. 

The facility constructed was a set of twin reverberation cham- 
bers made with plywood frame and ma<;onite interior. The cham- 
bers were shaped like a five-sided parallelepiped with an inclined 
roof to facilitate sound reflections in the interior. Opening 
was made at one of the side walls to f.t a composite panel 16" 
by 16" in size. The sound transmission area was 12" x 12". The 
chamber, having an internal volume of 48 cubic feet with corres- 
ponding surface area of 83 square feet was adequate for sound 
transmission loss measurement at input frequencies above 400 Hz. 

This part of the project has supported one graduate student 
to complete his MSME thesis titled "Design, Construction, and 
Performance of a Small Reverberation Chamber for Transmission Loss 
Measurements" (by C. Voorhees, 1974) and resulted in a journal 
publication co-authored by the Principal, Co-principal Investiga- 
tors and the graduate student (Applied Acoustics #9, 1975, pp. 
165-175). A copy of the reprint of the journal article is included 
in the Appendix. Preliminary result of the investigation was 
presented in the 87^*^ Meeting of the Acoustical Society of America, 
New York, NY, April 23-26, 1974. The abstract of the presentation 
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appears in the Journal of Acoustical Society of America Vol 55, 
Supplement, p. 513, 1974. 

The test facility was not put to use in the NASA Ames wind 
tunnels due to scheduling difficulties. 

4. Measurement of Transmission Loss of Composite Panels 

The sound insulation property of several composite material 
panels were measured. The panels were placed between the twin 
reverberation chambers. Bands of third octane sound were generated 
by speakers placed in the source room. The transmission loss for 
each panel was obtained by modifying the noise reduction across 
the panel, taking into account the absorptive power of the walls 
of the receiving room. The results for three composite panels are 
presented in Figure 1. These were graphite epoxy composite panels 
0.038" thick. Each panel was made of eight laminations with the 
graphite fibers orientated respectively at ±45", ±22.5° and 0-90° 
to one of the edges. Experimental results showed that the panels 
behaved very similarly to one with homogeneous structure such as 
comnon aluminum panels. In the upper audible range from 500 Hz 
to ultrasonic, the behavior of the panels is governed by the limp- 
wall law. A comparison of the transmission loss of a 0.039" graphite 
epoxy composite panel to one of 0.026" aluminum panel, both of same 
surface density are gi/en in Figure 2. Thus the composite panel, 
due to its high rigidity, may have an advantage in reduction of 
structural vibration, but its light weight have not resulted in an 
improvement of sound insulation characteristic. This finding should 
prove useful to NASA space vehicle designers responsible in reduction 
of cabin noise during space flight. 
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Experimental and Theoretical Investigation of Stress Wave 

Attenuation by Inclusions 

J. C S. Yang* and C. Y. Tsui+ 

University of Maryland, College Park, Md. 


The propagadoo oi an initially sharp cylindrical pressure pulse through a linear elastic medium with ladniions 
is analy^re^ both experimentally and analytically. In the experiment, tests were performed on plates with single 
and multiple circular inclusions embedded in a matrix of lower characteristic impedance. Sharp comprcssioa pulses 
were generated at an edge of the plate by two methods; the detonation of a charge of lead azide and the impalement 
of a sharp hujfe-edged ram impacted at the opposite end by a short projectile. Strain gages were mounted on 
various positrons of the plate to determine the attenuation of the transient stress. The qualitative analytical 
treatment b based on the methods of propagating stress £$cootinuities. Computer programs were written to 
numerically detenmne the changes in the shape of the leading wave front and the stresses immedbtely behind H. 
Experimental results for the attenuation of stress wave on steel-aluminum and steel-brass indusioo-iiiatrix 
composites compared very wdl with the computed analytical results. 


Introduction 

G rowing requirements for li^twei^t materials with 
great strength and stiffness, especially in the aerospace field 
on such projects as STOL, VTOL aircraft. Space Shuttle vehicles, 
heat shield, etc, have led to increased interest and research in 
the area of fiber composite materials. 

One of the most interesting properties of fiber-ceinforoed 
cc , >ite$ is their dispersive characteristics. In general, when a 
stress pulse is applied to liic boundary of a oomposile, the 
shape of the pulse is changed, and it **smears out** as it pro- 
pagates through the slab, with an attenuation of amplitude. 
Thus the composite appears to have greater dynamic than static 
strength. However, for dynamic loads, the simplest geometrical 
configuration may give rise to tensile stresses, even thou^ the 
loading is compressive. The load-carrying capacity of a body 
fabricated bn a fiber-reinforced composite is adversely affected 


Presented as Paper 72-394 at the AIAA/ASME/SAE 13th Structures, 
$ .. tural Dynamics, and Materials Conference, San Antonio, Texas, 
April 10-!2, 1972; submitted April 25, 1972; revision received 
S^tember 22. 1972. This work was jointly supported by NASA 
Ames Ktsearch Center under contract NGR-2I-M2-350; the Minta 
Martin Fund, College of Engineering and the Computer Science 
Center, University of Maryland. 

Index categories* Structural Dynamic Analysis; Structural Com- 
posite Materials including Coatings!. 

* Professor Member AlAA. 

^ Assistant Professor 


by bond failure at the interface bonds; in particular, it may be 
substantially smaller than the strengths of the materials con- 
stituting the composite, and relatively small tensile stresses may 
thus produce bond failure. In addition, fracture and spallation 
could occur due to the reflection at the back free surface of the 
slab, and internal stress concentration could occur due to re- 
fractions, reflections, and a focusing effect of the geometrical 
configuration. Few analytical and experimental investigations of 
such situations appear to have been carried out. In many cases, 
it is believed that the development of composites has proceeded 
with little appreciation for tne importance of microstnictural 
variables sudi as interface strength, moduius ratios, fiber size, 
etc, on the composite strength under stress-wave loading. Con- 
sequently, in order to have increased resistance to stress-wave 
fracture, a better understanding of the micromechanics of the 
fracture process must be developed. This approach is necessary 
to guide composite material development, and should result in 
large-scale practical benefits. 

In this paper, we are concerned with the propagation of an 
initially sharp cylindrical pressure pulse through a linear elastic 
medium with inclusions. Thv pattern of transient waves is 
extremely complicated, and for the analytical treatment, the 
attention is, therefore, focused on the changes in the shape of 
the leading wave front and the str^ immediately behind it. 

The basis for the present approach is the ray tracing and 
associated wave front analysts of geometrical optics. It can be 
shown that even in the case of dastic waves with two basic 
wave speeds in each material, the stress amplitude at the initial 
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wave front is given exactly by the law of geometrical optics 
associated with the propagation of irrotational waves. This 
follows from application of the theory of characteristics and is 
in accordance with the more common interpretation of geo- 
metrical optics as an approximation to the solution for oscillatory 
waves at high frequencies. This approach determines the stress 
tensor at the wave front, and can be extended also to determine 
stress gradients there. 

The propagation of stress discontinuities was discussed in 
considerable detail by Thomas. ‘ Much of the recent work in this 
area is patterned alter that due to Luneburg,^ for electro- 
magnetic discontinuities, and Friedlander,^ for the propagation 
of discontinuities in acoustic media. Propagation of stress dis^ 
continuities in elastic media, including the reflections and 
refractions of discontinuities at interfaces of different elastic 
solids, was investigated by Keller.^ An application of wave-front 
analysis to the problem of an inclusion embedded in an elastic 
medium was presented by Ting and Lee.^ Achenbadi, Hemann, 
and Ziegler^ presented an inter ^ting experiment in studying 
separation of the interface of a circular inclusion and the sur- 
rounding medium. Nariboli^ combined the ideas of the theory of 
singular surfaces with those of ray theory. Recently, Ben-Amoz* 
considered the problem of propagation of actual stress pulses 
through bound^ anisotropic media. 


Experiments 

The specimen used for the experiment consisted of 12 in.^ 
plates made of aluminum or brass with steel inclusions. The 
inclusions were shrunk fit in the plate with an interface pr^ure 
of approximately 200 psi. The steel-aluminum combination has 
a characteristic impedance ratio of 2.66 and a propagation 
velocity ratio of 1.03, and the steel and brass have ratios of 
1.27 and 1.39, respectively. Four different kinds of specimens 
were tested— i.e., one, two, three, and four cylindrical inclusions 
with the total area kept equal. A drawing of a single-inclusion 
specimen is shown in Fig. 1. In addition, experiments were also 
conducted on 12 in. x 12 in., 12 in. x 48 in., and 48 in. x48 in. 
aluminum plates without inclusion. 

The compressive pulse in the plate was generated by two 
methods: 1) The detonation of lead azide at the center of the 
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upper edge (The plates were placed in a vertical position and 
supported by rubber pads at the lower comers.); and 2) The 
impingement of a knife-edged ram (The specimens were sup- 
ported freely on rollers with the center of one edge resting 
against the ram. The other end of the ram was impacted by a 
projectile released from an air gun.) A schematic diagram of this 
experimental arrangement is given in Fig. 2. 

Both the lead azide charge and the impingement of the sharp 
ram produced cylindrical waves emanating from the edge at the 
point of contact into the plate. During the experiment, different 
amounts of lead azide and several impact velocities and projectile 
lengths were used. Pairs of strain gages were mounted on the 
opposite faces of the plate at various positions along the rays 
(see Fig. 1) to measure the transient stress. Each pair of gages 
are connected in the opposite arms of a wheatstone bridge 
circuit to cancel any bending effects. The output signal from the 
bridge circuit was fed into an oscilloscope. 


Analysis 

The wave fronts can be located by a direct calculation using 
Fermat's principle. To illustrate the procedure in computing rays 
and wave front, let us consider a ray AB in Fig. 3 which strikes 
the inclusion at angular position 0. For simplicity, distortion 
wave fronts are omitted and only the rays associated with a 
compresrive wave front are sketched. If a is the radius of the 
inclusion and b the distance from the applied pressure point to 
the inclusion, the coordinates of point B are given by 
xg = LasCOs 9 and yg = LAgrin^. where LAg is the length of the 
straight ray AB and is given by 

Lab - ib-^a) cos0-[(b+fl)^ cos^0-h^-2nh]*^^ 

Let r s 0 be the time at which the cylindrical wave initiates 
at X s 0. The time r required for the wave front to reach point B 
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Fit* 5 Refriciioo of two nciglibortet mys. 


is r *= where Ci is the speed of longitudinal waves in 

medium 1. ITie time r required for point C on the refracted ray 
BC to be reached by the wave fron* is 

f = Lab/c, -f Lg( /cj. or Lb( B c^If- 

where is length of the straight ray BC and is the 
speed of longitudinal waves in medium 2. The inciinatitm of the 
r^racted ray to the incident ray direction, 6. is given by Sneirs 
Law 

c,/c^ » sinC/$in0 (1) 

where C 0+atand ^ 

The coordinates of point C are given by .kc« v*+ 
LBcCos(d-»><$|and >’c = >*B'^LecSin(9+^). where the length 
can also be expres^ by = 2ncos^. 

A similar calculation determines the wave front after refraction 
out of the inclusion. The coordinates of point D are givoi by 

.Ko B .tc+Lco sin (0+24) and >d » >‘c+LcD^n(0+24) 
where JLcd* length of the straight ray CD, is determined by 
t = Lab/Cj +L ec/o**'^<iy^^i Leo = Cil— L ab— C|/c2Lbc* 

The wave fronts as expressed by the aforementioned equations 
are shown in Fig, 4 for various times. 

A complete analytical treatment of the stresses for the wave 
propagation problem shown in Fig. I is very difficult. However, 
the stresses at the instant that the wave fronts arrive can be 
obtained in a rather straightforward and simple manner by 
investigating propagating stress discontinuities. The analysis 
assumes linear elastic response of the materials of the inclusion 
and the matrix. In addition, it is assumed that the applied 
pressure instantaneously rises to a certain value and then rapidly 
decreases to zero. The stress discontinuity applied on the 
boundary generates surfaces of discontinuity whidi propagates 
into the interior of ihe continuous body in a manner which 
was discussed by other authors, most notably by Thomas* and 
Keller.^ Here we consider the problem in two spatial dimensions. 
When given the initial signal at a point on each ray, the 
variation of this signal along the ray is governed by 

p^R = const (2) 





Fig. 7 Repfodoclioa of strain bisiorics due to exploiavc loodug on a 
single-indiision steeLaluiramim spedmen, 0^ ray* 


Changes in intensity at the boundary between the matrix and 
inclusions are determined by the equivalent of Fresnel's formulas 
for elastic media.^ The pressure ratio at the interface can be 
written in the following form : 

p p(c,*-4c^^/3)c,* 3Mcosy 
p“pk,*-4c,»/3)c,* KM+LN 

where 

K =cosC^ + wsinC^^ j^I-^^sin (l-#isin^C)j 

\P J 

T. 

A? = |l — - — «sin*)siiiC-^nsiiiJcosJx 
\ P / Ci 

i-e-)T 


The quantities with a bar refer to the inclusion whereas the 
quantities without a bar refer to the matrix. 

Let and Pb he pressure at point B in Fig. 5 on the side of 
medium 2 and medium 1 respectively. 



Fig. 6 Pirasurt at wave froot for cylindrical indniofi* 
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The pre»urp can be deteonined by using Eq. (2) 
Pi-PbCfc/^Aj''* 

where is the p essure-at point 0 and b is the length of the 
straight ray AO, and the ndracted pressure ^ can then be 
determmed by u ng Eq. (3). Then^ Eq. (2) is used to relate 
points B and C 

Ai/Pa = 

where s the radius of curvature of the wave fronts at point B 
in mediun 2, and Ac radius of curvature of the vrave front 

at point C ;n medium 2. The radius of curvature, 
refracted wave fror ts at point B can be obtained by considering 
two neighboring points on the interface, B and B' as shown in 
Fig. 5. The angle of refraction at B' is 


W/dO(AC+A9) 

From Sneirsl .aw [Eq. (1)], we find d^/dC =» tan^/tanC. 
The followLtg relations are deduced from Fig. 5: 

.3C«BN/Lab. 

and 


(d^/dC)(AC+A^“ * FW/Kj 

Also BN = EB coSs. = Bff cos^ After substituting and 
simplifying, we obtain : 


^0 -cos^ 


pan 0 /cc^ + iV M 

[_tanf \Ub «/ 


I 


Snee the rays are $trai}^t lines normal to the wave front in 
both medium we obtain 


Ac ® 

The pressure pc can br determined from Eq. (3) where pc is the 
i>res5ureat point C afier *)eing refracted out of the inclusion to 
the matrix. Then using Eq (2) 

Fd/^c = 


Sen. 0.1 /vJD 


Sen. n. 05 v/an 

cm. 

I 

Sweep rate 10 Msec/cin 





The derivation of the radius of curvature Rc follows the same 
procedure as before: 


Rc^ cosi^^S) 


tan(^"f^) / cos^ _ ^ ll 

tan^ \ Rc a) aj 


-t 


Again, since the rays are strai^t lines normal to the wave front 
in each medium, Bp « 

Figure 6 shows an example of the pressure p as a function of 
position whoi the wavefront arrives. The r^ltsare shown in the 
form of contour lines of constant pressure. The material con- 
stants used are Cj/Ci « 1.39, and characteristic impedance 
ratio % 1.27. The pressure at point 0 is taken to be unity. 
Figure 6, when used in conjunction with Fig. 4, enables one to 
determine the pressure at the wave front for a given r. The 
wave fronts chained in Fig. 4 are alfo shown in Fig. 6, 
Oearly, one can continue this procedure to determine the wave 
fronts and the pr^ure at the wave fronts for more than one 
inclusion. 


Results 

The magnitude and shape of the experimental strain histories 
were plotted from tyfucal osdlloscope traces. A sample plot for 
the aluminum specimen with a sinj^e steel indusion is shown 
in Fig. 7. The result was due to the dc.onation of a charge of lead 
azide at point A (see Fig. I). The upper and lower traces 
correspond to the strain measured by the two pairs of strain 
gages which are located along the 0° ray, 4 in. and 10 ia from the 
point of detonation. The center of the 2-in.-diam inclusion is 
midway between the gages. 

A rise time of 4 psec was observed on both traces for the 
pressure pulse to reach its peaks G and H. The delay time of 
29 psec for the arrival of the pressure pulse between the gages 
correlated very well with the theoretically computed delay time 
of 29.5 /isec. The amplitude ratio of the peak strain, H/6, for 
the two strain gage measurements is 0.52. This compared very 
well with the amplitude ratio of 0.55 calculated from the ray 
theory. Another experiment was conducted on a specimen of the 
same dimension except without the inclus, ^n. The result is shown 
in Fig 8. The amplitude ratio of the peak strains. ///G. is 0.60 
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E = 3 

1 

0 

B53 

055 


aluminum 

F- 10 
E-3 

, 

0 

B42 

043 

11 

Ram 

II 

F-9 
E- 1.8 
F-9 

2 

0 

032 

0315 

impingement 

E-2.4 

3 

0 

0.34 

034 



F- 10 
E-2.4 
F- 10 

4 

0 

0.50 

056 


Aluminum 

E-4 
F- 10 

0 


B60 

063 



E-3 

0 



057 

III 


F-9 
E-24 
F- 10 

0 



0.49 



E- 18 
F-9.0 

0 
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as compared with a theoretically calculated value of 0i63. The 
delay time between the two gages again correlated well with the 
theoretically computed delay time. It should be noted that the 
amplitude ratio of the peak strains has been attenuated from 
0.60 to 0.50 for the plate with the inclusion. 

The strain variation at the location of each gage is also 
analyzed as a function of time. However, the calculation gets 
very complicated with the existence of the reflected and refracted 
waves. When there exists a large number of interfaces as in the 
cases of specimens with multiple inclusions or specimens with 
small overall dimensions, the calculation becomes virtually im- 
possible and some effective modulus theory must be used. The 
pulse reflection from the back free surface is also shown in 
Figs. 7 and 8 with the peak tensile amplitude /. The rather 
high peak amplitude J can be attributed to the multiple re- 
flections along the 0° ray and to the reflections from the sides 



Fig. 13 Strain histories tfur to eii|ilosivf loading on a two-hiclurion 
fteeLbrass spedmcn, 0° ray. 


of the square plate which arrived at the second gage approxi- 
mately, ^ fiscc after the initial detonation of the charge. Both 
peaks, I and J appeared during tests on plates of the same 
dimension with or without the inclusion ; and disappeared during 
tests on a similar plate but with a dimension of 48 in. x 48 in. 
as can be seen from Fig. 9. 

Additional sample results are presented in Figs. 10-14. 
Figure 10 gives the strain histories of the same steel-aluminum 
specimen along a 5° ray. Figures 11 and 12 are the results for 
a steel-brass plate with one inclusion along rays 0 « 0*" and 4^ 
Figure 13 presents the strain histories for a steel-brass plate with 
two inclusions. 

While the above are the results due to explosive loading, the 
strain histories of the steel-brass plate with a single inclusion 



Hg. 14 Strain kittorics 4m to impingement loading on a single- 
inclosion steel-brass specimen, 0" ray. 
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due to impingment of the ram against the plate at A (see Fig. 1 ) 
are presented in Fig. 14. The risetime has a much lai^ value 
of 10 fiscc compared to 4 iisec for the explosive loading. 

AJ^tional experiments using both methods were performed 
on plates with two, three, or four indusion& The results of all 
the experiments are summarized in Table 1. 

Table I is divided into three parts. Parts 1 and II refer to the 
two methods of pulse generation. Part 111 lists the theoretically 
computed peak ratio^ H/G for several selected gage locations. 
The results are based on Eq. (2) and are for homogoieous 
material The first four columns in Table 1 give, respectively, 
the composites tested, the location of the strain gages, the 
number of inclusions, and the identity of the rays. The last two 
columns give the average experimental results and the r^ults 
obtained by the ray theory presented in this paper. 


Conciiisiofi 

Experimental results on steel-aluminum and steel-brass 
inclusion-matrix composites using ex|^osives and ram impinge- 
ment compared very well with the analytically computed results 
using ray theory. 

With good agreement betwera the experimental and 
theoretically computed results, one can confidoitly proceed in 
applying the theory to fiber-reinforced composites to determine 
the optimum micro-structural variables such as modulus of 


elasticity ratio (fiber to matrix), relative fiber size, and orientation 
to increase resistance to stress fracture. 
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ISOLATION AND DAMPING 


1H£ MEASUREMENT OF DAMPING AND THE DETECTION OF 
DAMAGES IN STRUCTURES BY THE RANDOM 
DECREMENT TECHNIQUE 


J. C. S. Yang and D. N. Caldwell 
Mechanical Engineering Department 
University of Maryland 
College Park* Maryland 


A technique called random decrement has been developed, 
which makes possible the computation of damping values and 
the detection of damage In structures when only response data 
Is available. Damping ratios were computed using this technique 
for several modes of randomly excited panels, beams, and bones. 
These damping ratios compared satisfactorily to damping ratios 
which were computed from the power spectral density method. 
Standard randomdec s1gnatu<^e$ were established for all the 
structures. Damages were detected by bserving the changes In 
the established signatures. Notches f,n1ch simulated cracks were 
induced Into two of the beams. The effects of these notches on 
the beams' signatures are presented. 


INTRODUCTION 

When a physical structure Is subjected to 
random forces, Internal vibrations are set up 
within the structure as It absorbs the energy 
Imparted to It. In some cases, the energy 
Imparted may cause stresses which exceed the 
strength limitations of the structure and which 
may result in a failure cf the j»tructure. To 
prevent the failure of a structure, structural 
damping needs to be more thoroughly Investigat- 
ed. Small cracks and localized failures of the 
structural elements often have a significant 
effect on the vibration response characteristics 
of the structure long before they are signifi- 
cant enough to be visually detectable. Damping 
is one of the characteristics which changes* 
Damping is the means by which structures absorb 
energy and significant changes In damping are 
therefore very meaningful. The precise measure- 
ment of damping Is therefore highly important 
to the economic design and reliable analysis of 
large structures. 

Various attempts have been made to p«*ovide 
means for obtaining structural damping Inform- 
ation used In the design of structures and used 
In the monitoring of the response to the applied 
forces. Although various types of apparatus 
and techniques are available for measuring 
structural responses to randor^ vibrations, the 
data obtained Is usually so complicated that an 
observer cannot readily determine when a signi- 
ficant change In the structural response occurs. 
Most of these techniques are only suitable for 
use under controlled laboratory conditions and 
are of little use for structures In service. 


A simple, direct, and precise method Is need- 
ed for translating the structural response time 
history Into a form meaningful to the observer. 
The method chosen to examine the structure 
should be sensitive to changes In the natural 
frequencies an'^ the damping. Power spectral 
density has been considered, with damping meas- 
ured by the half -power point bandwidth method, 
but this was found to have a large measurement 
variance, especially when the bandwidth was 
small. In addition, when two modes are close, 
this method cannot be applied. Erroneous 
answers were obtained when assumed linear 
systems were actually nonlinear, a problem which 
could not be detected unless the Input was also 
measured. The autocorrelation function was In- 
vestigated as an alternative wherein damping 
data was obtained from the logarithmic decre- 
ment, The problem with the use of autocorre- 
lation signatures is that the level of the curve 
is dependent on the Intensity of the random In- 
put, and In a natural environment this can sel- 
dom be measured or controlled. If the structure 
is a linear system, the level changes can be 
compensated for by normalizing the curves, but 
If the structure is nonlinear (as Is often the 
case), a different signature will be obtained 
with each level of excitation. Therefore, cor- 
relation functions can only be used with linear 
systems and by knowing the Input. 

Another area that needs development 1$ the 
detection of crack Initiation and growth. Pre- 
sent methods of crack detection Include visual 
Inspection and acoustic emissions. Although 
acoustic emissions can detect flaws In assembly 
line comparisons. It Is highly unlikely that 
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under conditions of high ambient noise level 
such as Is encountered In aircraft flight that 
this method can be applied. In addition » It Is 
obvious that visual Inspections are useless 
when cracks develop within the Interior of a 
structural material. 

In the present report a technique called 
"Random Decrement" Is presented which makes 
possible the computation of damping values and 
the detection of damage In structures when only 
response data Is available. 

RANDOMDEC ANALYSIS 

A technique called "Random Decrement" anal- 
ysis has been developed and explored Initially 
by Cole Cl], which advances the state-of-the- 
art In the detection of crack development and 
extension In structures subjected to random 
excitation when only response data Is available. 
A brief summary of the method Is given here. 

The techniqe, or method of analysis, requires 
as Input the time history of the response of a 
system subjected to random excitation This 
time history is divided Into short segments, 
with each segment starting with an amplitude Y . 
This amplitude Is usually chosen to be some ^ 
fraction of the standard deviation of the sig- 
nal. The segments are ensemble averaged to 
yield the signature 6(t) which Is a time vari- 
ation of the parameter, say acceleration, 
strain, etc., composing the time history. The 
averaging process removes the effect of the 
random excitation and leaves a signature repre- 
senting the free vibration characteristics of 
the structure. The averaging process also 
yields a signature with an Initial slope of zero 
since the segments are divided equally Into 
segments of positive and negative Initial slopes. 
The method of extracting a randomdec signature 
Is shown In Figure 1. 

The degree to which the random excitation Is 
removed from the signature depends on the num- 
ber of segments averaged. Generally, good 
repeatability Is obtained with about 500 seg- 
ments. If high frequencies are being studied, 
say 5 KHz, a few seconds of data are needed to 
obtain an accurate signature. 

An li/iportant part of the analysis Is the 
filtering that Is done to Isolate certain fre- 
quency ranges. The time history Is usually 
bandpass filtered; so the resulting signature 
shows the characteristics of the structure for 
that frequency bandwidth. To obtain damping, 
only the natural mods of Interest Is analyzed. 
This effectively reduces the response to that 
of a single degree of freedom system. Once the 
signature Is obtained, damping can be found by 
the classical logarithmic decrement method. In 
addition, because of the fixed selection level, 
Yc, It was shown that the scale and form of 
this signature is always the same even though 
the random excitation may change intensity. 

This Independence of Input Intensity makes this 
method very attractive to use as a failure 


detector. 

To detect flaws, a standard randomdec signa- 
ture Is established for the undamaged struc- 
ture. Then, a flaw such as a fatigue crack 
will Introduce an additional degree of freedom 
for the structure and the frequency of the flaw 
mode will decrease as flaw size Increases. Thus 
when the flaw frequency couples with the modes 
within the analyzed bandwidth, the signature 
will change shape. An additional mechanism 
which causes distortion of the signature shape 
Is the frictional damping caused by the Inter- 
action of the crack faces. One Important 
characteristic of this method Is that It re- 
quires no knowledge of the structural excita- 
tion other than the fact that it If '^andom. 

EXPERIMENTAL PROCEDURE 

A. Analog Signal Recording 

Three different types of structures were 
tested: rectangular cross-section beams, bones, 
and thin square panels. Table 1 contains des- 
criptions of the specimens. Several different 
bones were tested. However, their first natu- 
ral frequencies were at about 2,000 Hz. The 
digitizing rate used was too slow to obtain 
randomdec signatures of sufficient smoothness 
to measure damping for this high frequency. 
Therefore, a mass was mounted on top of the 
bones to lower their first natural frequencies. 

Each beam was clamped horizontally In the 
center and mounted on an exciter. Acceleromet- 
ers were used to pick up the Input excitation 
and the output response of the beam to flexural 
vibration. An example of a mounted beam Is 
shown In Figure 2. 

Each bone was mounted vertically on a cylin- 
drical hollowed out aluminum base which was 
bolted to the shaker, see Figure 3. Accelero- 
meters were positioned inside the aluminum base 
and on the top end of the bone to pick up the 
Input excitation and to sense the response of 
the bones to the longitudinal vibrations. 

The panels were bolted with approximately 
clamped boundary conditions between two like 
reverberation chambers. The mounted panels 
can be seen in Figure 4. The panels were 
excited with a diffuse sound field with a sound 
pressure level of 105 dB re 0.0002 microbar. 

The excitation was provided by a loudspeaker 
located within one of the chambers. A conden- 
ser microphone was also placed within this 
chamber to record the signal which Is represen- 
tative of the random acoustical pressure vari- 
ations occurring on the side of the panel. A 
miniature accelerometer was used to sense the 
response. This small, light (0.4 gram) accel- 
erometer was utilized in order to minimize the 
effect of loading the structures. 

A random noise generator, connected to a 
power amplifier was used to drive the shaker 
and loudspeaker. Each excitation signal used 
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wa$ band-limited white noise. The white noise 
was low pass filtered at a cutoff frequency of 
about 3.500 Hz, This filtering was done to 
prevent the folding over of frequencies above 
the nyquist frequency during later digitization, 
and also to usefully utilize the tape recorder's 
full dynamic range. 

Before the Input and output signals were 
tape recorded, they were passed through Audio 
Frequency Spectrometers. These spectrometers 
were used to condition the signals and provide 
a rough look at the content of their frequency 
spectrums which were displayed by a Level Re- 
corder. These spectrums were used to identify 
resonances and examine the low pass filter 
characteristics. Subsequently, the Input and 
output signals were recorded on the tape record- 
er. 

6. Digital Signal Processing and Damping 

Measurement 

Because a digital computer cannot be used to 
store and operate on a continuous-time signal. 

It was necessary to change the analog signals 
Into discrete-time signals. The analog to dig- 
ital converter which was used had a sampling 
rate of 10,000 samples per second. 

Power spectral density of both the Input 
excitation and the output response was calcu- 
lated. Each power spectral density resulted 
from 18 averaged 1024 point fast Fourier trans- 
forms. Subsequently from these results the 
structure's transfer function was calculated 
and graphically displayed. From the transfer 
function plot, the natural frequencies and half 
power bandwidths of the resonances were measur- 
ed. The damping ratios were then calculated. 

The same data as above was also used to cal- 
culate the randomdec signatures. Before apply- 
ing the randomdec procedure to the data, the 
frequencies outside the narrowband time history 
of Interest were excluded. This was accomplish- 
ed by the convolution of a digital filter with- 
in the randomdec program with the data. After 
the narrowband data had been calculated, their 
mean and standard deviation were computed. The 
mean value was used to change the data values 
to have a zero mean. If this had not been done 
the signature would not have ended at zero as 
In Figure 1. A fraction of the standard dev- 
iation was used as the selected threshold level. 
Each time t'^e signal values were found to cross 
this level, a new segment was summed to pre- 
vious segments, thus calculating the randomdec 
s1v;iaturc. When the digital values of the sig- 
nature were computed, a plot routine within the 
program displayed the curve. From this decay- 
ing sinusoidal signature, the damping ratio for 
various numbers of cycles were calculated by 
the logarithmic decrement method. The damped 
natural frequency was calculated from an aver- 
age of the periods of oscillation. 


C. Detection of Damage 

To see the effect that damage might have on 
a structure's standard signature, notches which 
were assumed to simulate cracks were cut Into 
Beams 3 and 4 as shown In Figure 5. The 
notches were cut with a saw blade which was 
approximately 0.10 cm wide. Both beams were 
notched 1.27 cm from the supporting fixture. 

The notch In Beam 4 was cut Into the full 1 
Inch width of the beam to a depth of 0.05 cm. 
Signatures from the beam's response were com- 
puted with this notch and also for five more 
saw cuts which extended the notch to the total 
depth of 0.30 cm. 

Beam 3 was used to find the effect which a 
notch cut Into a beam's side woul:i have on It's 
signature. This was not expected to have quite 
the pronounced result on the signature as might 
be expected from Beam 4. However, torsional 
modes should be introduced by this type notch. 
Six damage signatures were taken for this beam 
from the effects of six 0.20 cm saw cuts, which 
made the ^ tal notch width 1.50 cm. 

RESULTS 

The computer generated power spectral den- 
sity and transfer function curves of the re- 
sponse and the excitation of a typical damped 
beam are showi^ In Figure 6 and 7. The narrow- 
band response ot one of the resonant peaks of 
the response was used to create the randomdec 
signature as shown In Figure 8. Comparisons 
of the transfer function and the randomdec sig- 
nature damping ratio results for all the speci- 
mens are given In Table 2. 

Figure 9 shows three randomdec signatures 
extracted from different sections of the same 
time history. Each of these signatures contain 
approximately 500 segments. These curves show 
the portion of the length of Beam 4's standard 
signature which can be considered stable enough 
from which to detect damage. Figures 10 and 11 
show the randomdec signaturts as beams 3 and 4 
are damaged. 

DISCUSSION AND CONCLUSIONS 

The limitations Imposed by the sampling rate 
of the digitizer and the digital filter charac- 
teristics at low frequencies considerably re- 
duced the number of damping ratic comparisons 
which were available. However, the results 
from Table 2 show that for most of the compari- 
sons which were made, the damping ratios taken 
from the randomdec signatures and transfer 
function plots differed by less than a factor 
of two. Furthermore, the few comparisons which 
did differ by more than a factor of two were 
associated with very low values of damping. 
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A comparison of values of the resonant fre* 
quencles obtained from the transfer function 
plots with those calculated from the periods of 
the randomdec sig. -tures show that all but one 
differed by le^ than 6 percent. Considering 
that the randomdec signatures were calculated 
from only response data while It was necessary 
to also use the excitation data for the calcu- 
lation of the transfer function, these results 
are very g,>od. This Is especially true for the 
panels since the excitation spectrum was very 
non-uniform. 

Figure 9 shows the results of extracting 
three randomdec signatures from different oor- 
tlons of the same narrowband time history. It 
can be seen that for approximately the first 8 
milliseconds all the signatures are almost 
Identical. This portion of the signature can 
be used as a standard from which damage can be 
detected. 

Figures 10 and 11 show the signatures result- 
ing from the damaged beams. Both figures shovf 
that the effects of the saw cut In each beam had 
a significant effect on the standard signatures. 

The satisfactory results obtained from the 
experimental comparisons Indicate that the ran- 
domdec technique Is a valid method of maasurlng 
damping. Also, the randomdec signatures of the 
damaged beams demonstrate its ability to detect 
the induced flaws. 
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Figure 2. Mounted Beam on Slu 



Figure 3. Mounted Bone on Shaker 



Figure 4. Mounted Panel In Chamber 
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Figure 5. Notched Beams 
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Figure 10. Damage Effects on Signature 
of Beam 3 Broadband Response 


Figure 11. Damage Effects on Signature 
of Beam 4 Broadband Response 
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TABLE 1. Descriptions of Beans and Panels 



Beams 


Number 

Naterlal 

Length 

(cn) 

Thickness Width 

(cn) (co) 

1 

Aluminum 

42.5 

0.124 2.54 

2 

Aluminum 

42.5 

0.325 2.54 

3 

Plexiglass 

42.5 

0.318 2.54 

4 

Plexiglass 

42.5 

0.432 2.54 

5 

Glass 

34.9 

0.635 2.54 

6 

Glass 

37.5 

0.635 2.54 


Panels* 


Number 

Material Thickness 

(cm) 

Lamina Lay Up 
(degrees) 

1 

Aluminum 

0.084 

Homogenous 

2 

Hercules 3501/ 
AS Graphite/ 
EpOJQf 

0.099 

0/90/0/90/90/0/90/0 

3 

Hercules 3501/ 
AS Graphite/ 
Epo)^ 

0.112 

45/-45/45/-45/-45/45/-45/45 

4 

Hercules 3501/ 
AS 

0.109 

22 . 5/-22 . 5/22. 5/-22. 5/-22. 5/ 
22.5/-22.5Z22.5 


*A11 panels are 30.5 cm x 

30.5 cm 


TABLE 2. Consolidation of Results 


Specimen 

Cut-Off 

Freqs. 

Transfer 

Function 

Randomdec 

Signature 

Comparison 

f 

i 


(Hz) 

f(Hz) 

C 

f^jHz) 

^Avg 


^Avg 

Beam 1 

350-600 

455 

.025 

45E 

.014 

.994 

1.79 

Beam 1 

750-1050 

901 

.013 

896 

.013 

1.004 

1.00 

Beam 2 

250-550 

350 

40 

CM 

o 

387 

.003 

.982 

8.67 

Beam 2 

800-1300 

1060 

.008 

1076 

.001 

.985 

8.00 

Beam 3 

300-650 

433 

.032 

441 

.016 

.982 

2.00 

Beam 3 

650-1100 

868 

.033 

879 

.022 

.987 

1.50 

Beam 4 

400-900 

613 

.022 

622 

.029 

.985 

.76 

Beam 4 

. 900-1500 

1225 

.017 

1225 

.020 

1.000 

.85 

Beam 5 

900-1300 

1065 

.019 

1099 

.001 

.969 

19.00 

Beam 6 

800-1100 

868 

.004 

909 

.002 

.955 

2.00 

Bone 1 

500-1500 

985 

.019 

1039 

.019 

.948 

1.00 

Panel 1 

135-190 

166 

.051 

164 

.036 

1.012 

1.42 

Panel 2 

140-190 

175 

.017 

173 

.019 

1.011 

.89 

Panel 3 

190-250 

224 

.040 

193 

.042 

1.160 

.95 

Panel 4 

190-270 

215 

.061 

200 

.039 

1.0/5 

1.56 
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RBSPC^SE OF PLATES TO RANDOM PRESSURE LOADING 
FROM A TURBULENT BOUNDARY LAYER 


William T. Messick* 
Jackson C. S» Yang** 


The study of structural response to turbulent boundary layers is 
very important for determining fatigue life and noise transmission to 
the interior of aerospace vehicles. To gain insight into the problemr 
the response of clamped plates to a cross spectral density function 
representing the pressure loading from an attached turbulent flow is 
studied. As a preliminary study to determining the response of composite 
material plates, an isotropic plate is modeled with finite elements and 
its response is ccmipared to theoretical and experimental results. The 
modal frequency and ramdom response ligid format number 11 of KASTRAN was 
used for determining the natural frequencies, nodal patterns, and 
displacement power spectral densities. 

The plate that was modeled has clamped edges, is made of steel of 
density 0.27 pounds per cubic inch and measures 4.0 inches in the flow 
direction by 2.75 inches vide by 0.015 inches thick. A grid was chosen 
that has nodes at the quarter and mid points of the plate for con^rison 
with references 1 and 2. To keep computer running time to a minimum, a 
relatively coarse 6x8 grid was chosen. The CQUAD2 element with the 
coupled mass (COUPMAS5) option was used for the finite element model. 
Constraining the membrane decrees of freedom, the rotations about the 
normal, and the clamped edges yields a problem with 147 decrees of freedom. 
The eigenvalues in the frequency* range of 0 to 3000 cps were ceterminec 
using rigid format 3 in NASTRAK and compared with these calculated from: 



where v is Poisson's ratio 
£ is the modulus 
p is the density per unit area 
h is the plate thickness 
a is the plate width 

th 

and X. is the characteristic value of the i mode and a fxinctxon 
^ of the plate aspect ratio. 


* Mechanical Engineer, Naval Ordnance Laboratory, White Oak, Md. 

♦♦Professor, Mechanical Engineering Department, University of Maryland 



The computed NASTRAN values and the error associated with each mode based 
on the characteristic values tabulated in reference 3 are shown in Table 1 
The agreement is very good for the coarse mesh. 


Mode Number 
(ro-n) 

NASTRAN Natural 
Frequencies for 
a Clamped 8xB 
Grid (cps) 

Error 

1-1 

542.9 

0 . 1 % 

2-1 

856. U 

.03% 

1-2 

1330.0 

x.i\ 

3-1 

1404.2 

i.?% 

2-2 

1611.5 

0.3% 

3-2 

2087.2 

-0.7% 

4-1 

2210.4 

5.0% 

1-3 

2551.5 

2.3% 

2-3 

2791.3 

0.7% 

4-2 

2798.2 

1.5% 


Table 1 - Natural Frequencies Obtained From 

NASTRAN 8x8 Grid for a Clamped Plate 


The pressure fluctuations that occur on a flat surface due to a 
turbulent boundary layer may be expressed in terms of random variables 
which are spatially and teiriporarily correlated. For homogeneous 
turbulence, the cross spectral density depends only on the separation 
distance of two points. Furthermore, the cross correlation coefficient 
is separable in terms of streamwise and cross flow directions so that the 
cross spectral density is expressed in the form: 


Sp(o)) 



( 2 ) 


where and C 2 are the streamwise and cross flow separation, respectively, 
of two points, u is the circular frequency, and t)_ is the narrow band 
convection velocity. |p (C^,0,u) | and |p ( 0 ,C 2 >(*>}{ are correlations 



coefficients for the streamwise and cross-flow directions, respectively, 
which depend on the type of flow* This is the general form of the 
fxinction that was input to NASTRAN using the RANDPS bulk data card. 

The general form of the RANDPS input is 


g(w) (X + iY) 


( 3 ) 


where X and Y are constants* However, for the function of interest, 

X and Y should be functions of frequency. As a result, a separate set 
of RA17DP5 cards had to be generated for each frequency and a separate 
computer run made for each one* 

For the subsonic flow case, the cross correlation coefficients have 
been shown by Bull (reference 4) to be 


P^(Ct fO,w) 
P 1 


, g-O.l 

= e--037|Ci|/6* 


if u6*/U > 0.37 

c 

if idSVu < 0.37 
c 


(4) 


0p(C,C2.«^) = if U2I 2 log + 5.4 


45] 


* C.28 + 0.72 e 


-0.547U2I/6*- if 


ICg! < c 


where D = (0.59 + 0.3 e 
c 

U is the free stream velocity 


atnd 6* is the displacement boundary' thicteess 


C 

(5) 


Using the boundary layer thickness of 0.179 inches for the Mach number 
0.3 case (reference 4) the cross spectral density values were evaluated 
for 9 frequencies in the 0-3000 cps range. Only the modes in this 
frequency range were used for the modal solution. A structural damping 
coefficient of 0.009 was used in order to compare the results with the 
modal analysis in reference 1. The KASTRAN results, the modal analysis 
results of reference 1, and the experimental results of reference 2 are 
shown plotted in Ficures 1 and 2 for the displacement power spectral 
density of the panel quarter point and mid point, respectively. The 
running time on the University of Maryland UNIVAC 1108 using KASTRAN 
level 15.1.0 was 463 CPU seconds for the eigenvalue analysis and the 
displacement power spectral density at one frequency. With a checkpointed 
tape of the eigenvalue analysis, each succeeding run time was 149 seconds 
per frequency. 



The agreement of the NASTKAN results with the modal analysis of 
reference 1 is excellent considering the relatively coarse mesh. There 
is a discrepancy at the panel midpoint for the ninth natural frequency 
but the difference is less than an order of magnitude. The agreement 
with the experimental values is also very good. It appears that with an 
additional number of NASTKAN calculations, the root mean square response 
of the panel can be adequately predicted. However, although the displace* 
ment power spectral density at a particular frequency may be obtained at 
a reasonable cost, the generation of a sufficient number of points for 
the root mean square displacement becomes an expensive proposition. To 
help decrease the running time it would be advantageous to allow the real 
and imaginary parts to be separate functions of frequency and let these 
functions be evaluated at a ntmber of frequencies for one computer run. 

Future studies will investigate the effect of filament orientation 
on the random response of composite plates. An orthotropic constitutive 
Tiatrix will be input and a finer mesh used for the NASTKAN analysis. 
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Figure 1 - Displacement PSD Response at Panel Quarter Point 
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THE DESIGN OF SMALL REVERBERATION CHAMBERS 
FOR TRANSMISSION LOSS MEASUREMENT 


Chung Y. Tsul Carl R. Voorhees* and Jackson C. S. Yang 

Department of Mechanical Engineering, University of Maryland, 
College Park, Maryland (£/S^) 


SUMMARY 

This paper describes the design^ construction and performance evaluation of a small 
test jacilityfor transmission loss measurement. The low-cost test facility — in the form 
of a twin pentagonal parallelepiped — is made of laminates of masonite^fibreglass^ and 
plywood and is particularly adaptable for transmission loss tests of small panels. 



INTRODUCTION 

The accurate assessment of the sound insulation properties of panels and partitions 
is an important area in acoustics, as is indicated by past research. This aspect of 
research took on a new importance recently with the revelation that a quiet 
environment is imperative for healthy living and efficient working. In order to protect 
privacy, as well as prevent hearing damage to the public, legislation and ordinances 
against excessive noise exposure have been promoted by all levels of government. 

Sound insulation data are essential to architects and noise control specialists in 
their planning for an environment with comfortable noise levels. Economic 
considerations dictate that this desirable environment be designed into new 
buildings and structures, rather than be treated later as an "add on' feature. These 
data are also required when remedial modifications are made to existing designs. 

The sound insulation property is most commonly obtained by the two-room 
method, especially for panels and partitions. In this method, the test specimen is 
placed in the common wall between the source room and the receiving room, and the 
noise reduction between the rooms is recorded experimentally. The result is 
expressed by the transmission loss, TL, defined by: 

TL«10log,o(l/T) (1) 

* Present address; Applied Acoustics Section, National Bureau of Standards, Washington, DC 20234 
(USA). 
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where x is the ratio of the transmitted energy to the incident energy. Measurement 
procedures and precautions can be found in recommended standards such as ASTM 
£90,’ published by the American Society for Testing and Materials. 

According to the standard, the transmission loss is given by: 

TL = NR + 10Iog,oS- 10!og,ov42 (2) 

where NR is the noise reduction between the source room and the receiving room, 5 
is the area of the sound transmitting surface of the test specimen and A 2 is the total 
absorption cf the receiving room expressed in units consistent with 5. 

The derivation of this equation :s based on many assumptions. The most 
important is that the sound field in the rooms be sufficiently diffused. This is not 
surprising as, in the design of rooms for measurement, the attainment of a diffused 
sound field has caused the most concern. The design of small rooms is no exception, 
as is reported in this paper. 

Many methods have been devised in order to achieve a relatively high degree of 
diffusivity. Some of the more important are outlined below. 

Rotating reflector 

A motoriswd rotating reflector placed in a room effectively creates a variable 
boundary of the room which results in oscillation and the appearance and 
disappearance of modal patterns for a constant sound source. The circulation of the 
air in the room also breaks up the modal patterns to achieve a better diffused sound 
field. 

Reflecting surfaces 

Highly reflective surfaces are added to the rooms by hanging objects from the 
ceiling. This increases the sound reflecting area for a given room volume without 
appreciably increasing the room absorption, thus enhancing sound diffusion 
through multiple reflections. 

Non*para!le! waiis 

This enables the sound waves in a room to impinge against the walls more 
randomly at various angles of incidence. Given a constant volume, the total number 
of normal modes for an irregularly shaped room is very close to that for a room with 
parallel w' ' s. However, the degeneracy of modes due to the overlap of modal 
patterns in an irregularly shaped room is much less and, consequently, a better 
diffused sound field is obtained. 

Multiple frequency sound source 

A source such as a warble tone or a third-octave band noise has energy over a wider 
frequency range than a pure tone. This source is capable of generating a series of 
modal patterns, each pertaining to a narrowly grouped frequency bnnd. The 
applicatioii of such a source in the measurement of transmission loss is an 
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operational practice and is actually required in the .A6fM E90 procedures. It is. 
however, unrelated to the design and construction of the chambers. 

There exist many rooms^'^ which have the desirable sound diffusion 
characteristic. These rooms are very reverberant and are commonly referred to as 
reverberation ch imbers. Located at various acoustical laboratories associated with 
research organisations, universities and manufacturers of acoustical materials, these 
rooms are also used for the measurement of power outputs of sound sources and 
sound absorption coefficients. For the m .asurement of transmission loss, however, 
two adjacent reverberation chambers are required. 

Due to the fact that acoustical characteristics, including transmission loss at 
frequencies as low as 100 Hz and under^ are required for full-sized building 
components such as doors and wall assemblies, these rooms are generally large, 
ranging in size from 1000 to over 10,000 ft^. The walls are generally made of massive 
concrete, I to 3 ft thick. There are inherent drawbacks with these large rooms. The 
first of these drawbacks is one of economics. The building and equipping of the 
rooms is costly and generally regarded as a major investment for a good acoustical 
laboratory. The second and more critical drawback is the acoustical permanence of 
the rooms. Because of their large size, the acoustical characteristics of the rooms are 
difficult to change. Therefore, for a given test facility, it is nearly impossible to assess 
the effects of the many important parameters which would influence the results of the 
measurement. In the same way, it would serve no useful purpose to reconcile for 
measurement discrepancies results from non-identical test facilities, even though the 
method of measurement anc' :he specimen were identical. A large test facility is 
generally unfit for testing specimens which are small compared with full-sized 
building partitions or doors. This inadequacy is particularly significant for 
transmission loss tests as the sound energy must be transmitted from the source 
room to the receiving room through a tunnel formed by the thick walls. The 
specimen will not be exposed to sound through a wide range of angles of incidence. 
Errors in measurement will also be larger as the size of the specimen becomes smaller 
due to the increasing proportion of transmitted energy through the common wall. 

Small test chambers are justified because they avoid the above drawbacks. It is 
possible to quickly and inexpensively change any room parameter in order to 
investigate its influence on test results. By eliminating the tunnel effect, these 
chambers are ^iso more adaptable for testing small panels and components. The 
latter advantage becomes even more attractive in aerospace applications as of 
the panels and components are inherently small. In some cases, it is required to test 
new materials for their acoustical properties when only a small sample of the 
material is available. The drawback of small test chambers lies mainly in their high 
cut-off frequency, which means that it is not possible to obtain reliable sound 
insulation information in the low frequency range. Whether this constitutes an 
unacceptable handicap is dependent on the intended application and must be 
resolved before the measurement is undertaken. 
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DESIGN AND CONSTRUCTION OF TEST FACILITY 

The test facility to be constructed was designed to yield reliable transmission loss 
data for small panels. The internal dimension of the lOoms was limited to between ? 
and 4 ft at an edge. It was decided, because of the size, that it was not feasible to have 
motorised rotating reflectors, and that a pair oi twin chambers made of plywood and 
masonite laminated boards would be the easiest type of facility to construcL To 



increase sound diffusivity. non-parallel vertical walls in the form of a pentagon with 
an inclined roof and a level floor were adopted. It can readily be appreciated that in 
choosing this particularly simple configuration, correct dimensions of the wall 
boards, especially at the edges, can accurately be obtained through geometrical 
developmenL thus ensuring good acoustical seals. The shape of the twin parallele- 
piped facility, which is shown in Fig. I, very closely resembled those at the Kobayasi 
Institute in Japan.^ A schematic view of the test chambers' is given in Fig. 2. 
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The common wall and the floor of the rooms were made of 1 -7S in laminate of 
two 0- 7S in plywood and one 0-2S in masonite boards. The other walls am* :he ceiling 
of the rooms were made of 1 in laminate of plywood and masonite, reinforced by a 
wooden gnd of one-by-three ribs. A 0-25 in outer masonite sheet enclosed a 3 in 
lining of Owen-Coming 703 Fiberglas insulation. The laminated boards were 
assembled and all joints were glued and fastened with screws. Acoustical sealant was 
applied along all edges to ensure a ught seal against noise leakage. 

The 22 in^ x 4-Sin frame for mounting the test panel was constructed of 
hardwood pieces of 3 in x 4-S in cross-section, was rigidly attached to the receiving 
room and protruded through an opening in the common wall (Fig. 2). A panel size of 
I6in X I6in could be fitted into the frame and clamped securely between 2in steel 
facing plates, leaving an exposed area of 12in x 12ia 



Fig. 2. Schematic vievv of the twin reverberation chambers. 



To facilitate access to the test panel and the measuring equipment a 28 in x 34 in 
opening was made in a rear wall of each room (Fig. 2). The opening was closed by a 
door of the same construction as the wall and was sealed on the inner and outer edges 
with a 2 in aluminium flange and rubber gasket. Electrical wirirg was fed through a 
side wall at a location near the ceiling and the common wall (Fig. 2). The small hole in 
t!ie wall was subsequently plugged with cork and acoustical sealant 
A microphone traversing mechanism was designed to carry a microphone 
continuously over a linear path. When activated the unit traversed at approximately 
8 in/min and automatically reversed direction at the end of each crossing. The 
tra ersing guide could be mounted at two difietent ceiling positions in each room 


ORIGINAL PAGE iS 
OP POOR QUALITY 




170 


CHUNG Y. TSUL CARL R. VOORHEES. JACKSON C S. YANG 


with the microphone suspended at two heights to provide four different microphone 
paths (Fig. 3). 

In connecting the twin chambers, the mounting frame was pushed into a similar 
opening in the common wall of the source room and pressed against a resilient stop 
at the inner side. The common spacing between the rooms was about I in and was 
filled with Owen-Coming 703 Fiberglas. The small gaps between the common walls 
and the mounting frame were also filled with acoustical sealant This arrangement 
allowed no rigid connections, thus discouraging flanking and structural 
transmission between the two rooms. 



Fig. 3. Microphone traverse paths and measurement stations. 

The completed test facility consisted of a pair of twin chambers each having a 
volume of 48ft^ and a surface area of 83 ft^. 


TEST FAOUTY PERFORMANCE EVALUATION 

Several tests were conducted to evaluate the performance of the facility. These tests 
established the lower frequency limit and the minimum sound pressure levels in the 
receiving room in order to yield valid transmission loss data for the specimens tested. 

Diffusivity measurement 

At various third*octave excitations, sound pressure levels along microphone 
traverses were recorded in order to examine the spatial variations. Two 
representative sets of results over four different traversing paths are shown in Fig. 4. 
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These results were obtained at centre frequencies of 31SHz and 400Hz. Sound 
pressure level variations along the same patte at other centre frequencies were found 
to be less significant for excitations higho- than 400 Hz and were more significant for 
excitations with centre frequencies lower than 3 1 5 Hz. These resulu indicate that the 
lower cut-off frequency for the test facility is at about 400 Hz. This value is slightly 
lower than the derived 500 H4 according to an equation in ASTM E9Q, to assure the 
achievement of an adequate number of room modes. 


J- 5' H 



315 Hz' 


H 3* H 



400 Hz 


Fig. 4. Sound pressure level vanations at 315 Hz and 400Hz. 


Reverberation time measurement 

The reverberation times of each room were determined from the decay curves(Fig. 
5) of third-octave excitations ranging from 400 Hz to 4000 Hz. ASTM C423 
requirements^ were followed. In addition, the absorption coefficients of the walls of 
the rooms were obtained from the Sabine equation.’ The calculated results arc 
shown in Table 1. Each number in the Table is the average of 20 micn^hone positions 
shown in Fig. 3. In most cases {see Table 1) the absorption coefficients meet the 
ASTM E90 recommendation of 0 06 for good sound reflections. The constant decay 
slope of the curves in Fig 5 which were obtained at diflerent stations also indicates 
that a good diffused sound field exists at 400 Hz. 
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Fig. 5. Decay curves at 400 Hz. 



Field transmission loss measurement 

Field transmission loss measurement provides the noise reduction between two 
adjacent rooms separated by a common wall This concept was applied to the present 
case in the establishment of the minimum sound pressure levels in the receiving rcom 
o\^er the background noise levels surrounding the test facility for reliable 
transmission loss measurements. In this case, the two adjacent rooms were the 
receiving room and the laboratory which housed the test facility; the common wall 


TABLE 1 

REVERBERATION TIMES AND ABSORPTION COEFRaENTS OF THE TEST FAnUTY 


Third-^KUM 
centre freq, {Hz) 

400 

SOO 

630 

800 

1000 

1250 

1600 

2000 

2500 

3150 

4000 

T»o Source rooms 

04S0 

(1550 

0510 

0470 

0450 

0400 

0 470 

0450 

0460 

0430 

0460 

T^ Receiving rooms 
a Source room 

0430 

0484 

0503 

0503 

0 482 

0486 

0481 

0448 

0419 

0404 

0408 

0061 

0050 

0054 

0 058 

0061 

0068 

0*058 

0*061 

0059 

0*063 

0^59 

a Receiving room 

0064 

0056 

0054 

0 054 

0^057 

0 056 

0 057 

0061 

0 065 

0068 

0K)67 
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was the set of walls* ceiling* and floor of the receiving room. In the measurement 
ASTM E336 requirements® were followed and the results are given in Table 2. 
Therefore, an ambient noise level of 40 dB in the laboratory will result in a 
background noise level of 3dB in the receiving room; for reliable transmission loss 
data, the recorded sound pressure level should be at least 10 dB higher or above 
13dB in the receiving room. 


TABLE 2 

FIELD IHANSMISSION LOSS OF THE RHrEIVING ROOM WALLS 


Octave centre freq. {H:) 

250 


1000 



IJUboratory 
SpL Laboratory, dB 

018 

80 

025 

») 

(V30 

80 

027 

75 

026 

75 

iMkeceivmg room 
SpL Receiving room, dB 

0052 

(M152 

(H)63 

0063 

C062 

58 

57 

57 

47 

43 

Field transmission loss, dB 

54 

35 

34 

39 

43 



Flanking transmission loss measurement 
This is a measurement of the transmission loss of a panel of very high sound 
insulation resistance. The result indicates the degree of flanking transmission 
through paths other than the test panel. For a panel under actual tesL the 
transmission loss measurement should then be somewhat less in order to ensure that 
most of the sound energy is transmitted through the panel. The flanking 
transmission measurement then established an upper limit for valid transmission 


50- 



Fig 6 Banking transmission loss of test facility and region of reliable measurement results. 
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measurement results for the test facility. In this test, a 0*5 in steel panel was used and 
the results are given in Fig. 6. 

Based on the evaluation tests, a region of reliable transmission loss measurements 
for the test facility can be obtained. This region is bounded by hatched lines (Fig. 6) 
on the left at 400 Hz and from above by a line approximately parallel to the flanking 
transmission loss curve but set arbitrarily lOdB below it. Transmission loss 
measurements for a panel within this region are valid and reliable. 



The test facility was used to obtain transmission loss data for two 16in x 16 in 
panels with the results shown in Fig. 7. Result 1 refers to a 0-033 in aluminium panel 
of surface density of 0*46 Ib/ft^. This curve is compared to a field incidence curve 
defined by:’ 

(TLU^-(TL)o-SdB (3) 

where (TL)o is the transmission loss calculated from the normal incidence limp-wall 
law. The result exhibits the well-known mass-controlled behaviour in the frequency 
range of 400 to 10,000 Hz. Result II refers to a 0*039 in graphite epoxy composite 
panel with 45 ° fibre orientation and surface density of 0-36 Ib/ft^. The panel exhibits 
stiffness-controlled effects in the frequency range 400 to 1000 Hz due to its high 
modulus of elasticity but reverts to mass-controlled behaviour in the frequency 
range 10(K) to 10,000 Hz. Due to the smaller surface density, the transmission loss for 
the composite panel is some 2 dB lower than that for the aluminium panel. If the 
thickness of the aluminium panel were decreased to 0*026 in, the surface density 
would be comparable with that of the composite panel and the transmission loss for 
the two panels would virtually be equal. 
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In summary, it is possible to conclude, based on results obtained, that the 0*039 in 
composite panel tested is as effective a sound barrier as a 0*026 in aluminium panel in 
the frequency range of approximately 1000 Hz to 10,000 Hz. When subjected to a 
random acoustical loading, a composite panel appears to react similarly to those 
made of more conventional materials, notably in the mass-controlled frequency 
range. Because of this, the advantage of the light weight of the material is offset by 
higher noise transmission. 
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